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Abstract 

Wave chaos is demonstrated by studying a wave propagation in a periodically corrugated wave- 
guide. In the limit of a short wave approximation (SWA) the underlying description is related 
to the chaotic ray dynamics. In this case the control parameter of the problem is characterized 
by the corrugation amplitude and the SWA parameter. The considered model is fairly suitable 
and tractable for the analytical analysis of a wave localization length. The number of eigenmodes 
characterized the width of the localized wave packet is estimated analytically. 
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INTRODUCTION 



Wave chaos phenomena reveal themselves in underwater acoustics. It is known that 
a wave propagation with the wavelength A in a long range-dependent wave-guide can be 
described under certain conditions by the parabolic equation in the limit of a small-angle 



propagation 



2j. This equation corresponds formally to the Schrodinger equation with 
an effective semiclassical parameter being of the order of X/L, where L is the characteristic 
size of inhomogeneities. It has been shown that this issue of quantum chaos is relevant to 
the phenomenon of long-range acoustic wave propagation in the ocean, where a ray-based 



approach is widely used in schemes of ocean acoustic monitoring 



numerical evidence 



a, |6j|. There is 



3, la, |9| that the phenomenon of ray chaos plays an important role in 
long-range sound transition. Wave propagation in the inhomogeneous wave-guide media 
can be considered from the point of view of nonlinear Hamiltonian dynamics, where the 



loj. One of 



consideration can be reduced to a corresponding problem of ray dynamics 
the most interesting features of ray dynamics is the possibility of the dynamical chaos of 
rays when a wave-guide carries longitudinalperiodic inhomogeneities. Unlike ray dynamics, 
where wave features are omitted, the important specific features of wave propagation can 
be taken into account by considering the parabolic equation for a small-angle-propagation 



111 ] . The relation between the wave equation and the corresponding ray dynamics was 



the issue of extensive studies 
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In this paper, we consider a model of the wave propagation in a range-dependent wave- 
guide. When X/L —>■ 0, i.e. in the short wave approximation (SWA), the underlying 
description can be reduced to the ray dynamics which can be chaotic due to range-dependent 



perturbations. The nature of this chaos is similar to the chaotic billiards 
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The formulation of the problem for the wave transmission in the framework of geometrical 



optics 



151 ] is relevant to the ocean monitoring problem where the 



parameter X/L is small: 



6) . However, a principal 



X/L ~ 10~ 3 - 10~ 2 , and an accuracy of the SWA seems fairly high 
difference between wave-like and ray-like descriptions, which may be negligible for the stable 
ray dynamics, has a strong impact in the case of chaotic rays. A manifestation of ray 
chaos in waves dynamics is similar to the effects of quantum chaos. We can indicate at 
least two such important manifestations: a breaking length x\ along the range-dependent 
(longitudinal) direction x, and a localization length which characterizes the cut-off of the 
number of eigenmodes An a of the propagated wave packet along x. The first one, x\, is 



16] , sometimes called the Ehrenfest 



similar to the breaking time in quantum chaos 
time. It shows an accumulation of small wave effects due to the chaos of rays, which leads 
to the breaking of the SWA for x > x\. The second manifestation of chaos related to 
localization of a wave packet in the momentum space and the corresponding cut-off of the 
number of eginmodes An\. This second manifestation of chaos is dynamical localization 
similar to the Anderson localization in solid state. 

The goal of this paper is to demonstrate the second manifestation of wave chaos using 
a fairly convenient model of wave/ray propagation in a periodically corrugated wave-guide 



(Fig. 1). The control parameter of the problem can be written as a product 



e/L = (2e/A) • (A/2L) 



(i) 



where the first multiplier characterizes the dimensionless amplitude of the perturbation 
(corrugation), while the second one is the SWA parameter. The problem, presented in the 
paper, is tractable for the analytical analysis, and it makes possible the calculate of the 
number of the localized wave modes An x for the two limiting cases both A ^> e and e A. 
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THE SEMICLASSICAL PARAMETER 



We first define the semiclassical parameter. Let us specify the general properties of 
solutions of a 2D wave equation in a homogeneous media bounded by two reflecting surfaces 
y = and y = yo + £<f>(x). Here 4>{x) is a corrugation function with the period 2%/ko, 
and the amplitude of corrugation is e, while yo = L is an average transverse size (or depth) 
of the channel (see Fig. 1). In what follows, we also consider that a small parameter is 
s/yo <C 1. The notation y Q is chosen to stress that the characteristic size has the transversal 
direction. The solution of the wave equation can be obtained in the following way. Assume 
that the time dependence of a wave field is harmonic, 



Applying conformal transformation to map the domain x G (—00,00), y G [0,yo + £ 4 > ( x )} 
to the domain x' G (—00, 00), y' G [0,yo] and omitting primes for shortness, one can write 
down in the first order of e/yo the stationary wave equation in the form: 



where A is the 2D Laplace operator, and u is the frequency of a sound wave with the velocity 
c s . It is obvious that the coefficient A(x, y) is a periodic function of x with the period of the 
corrugation function. Thus, due to the Floquet theorem, the solution is cast in the form 



where B(x, y) satisfies to the periodic B(x+27r/k , y) = B(x, y) and, e.g., Dirichlet boundary 
conditions, B(x, 0) = B(x, yo) = 0. In the zero order approximation, we obtain from ([2]) the 
following dispersion relation 



ty(x,y,t) = exp(—iut)^(x,y). 




(2) 



*(x,y) = e tkx B(x,y), 



(3) 



u 2 =u 2 m {k) =c 2 s [k 2 n +p 2 m }. 



(4) 
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Here p 2 m ~ (ir /yo) 2 m 2 is the discrete spectrum of transverse wave numbers (as a result of 
the boundary conditions) with the corresponding eigenf unctions B(x,y) = b m (x,y), where 



m = 0, 1, . . . 18]. The dispersion relation for the fixed u is a constraint for both the 
longitudinal k m and the transversal p m wave numbers. Therefore the solution is a finite 
sum on m in a range < m < m max : 

*(x,y) = e ikmX b m (x,y) . (5) 

m=0 

Since p m ~ ixm/yo for any type of the boundary conditions and the maximal momentum is 
Pmax = w/ca, we obtain that 

m max ~ ujy /7rc s = 2y /X. (6) 

Hence the validity of the geometrical optics limit (or ray dynamics) is determined by the 
condition m max = 2y Q /X ^> 1. In this approximation, an angle between a ray and the 
longitudinal coordinate is 

a m = arctan(p m //c m ). 

In the limit of small angles and using that k m w ou/c s , it gives the following SWA condition 

a m ~ i — = m = hm, (7) 

Km ^2/0 

where the inverse number of half-waves in the transverse direction of the wave-guide 

h = X/2y (8) 

is the semiclassical parameter. Transition to the classical or ray dynamics is therefore the 
following limits h — > and m — > oo. 
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THE ULAM MAP 



In what follows, we present a simple example of a ray dynamics quantization, where the 
non- vanishing values of h play an important role. We consider conditions where the limit of 
small angle propagation is chaotic. The ray propagation is characterized by two variables, 
namely, the dimensionless longitudinal coordinate x and the angle (or dimensionless mo- 
mentum) a. We suppose that inside the channel a ray undergoes complete refractions from 
the boundaries. Therefore the relation between the variables (a,x) for any two consequent 
bouncings, for instance, n and n + 1 (say, from the bottom boundary) are determined by a 
geometrical consideration which is presented in Fig. la,b. These relations form the following 
map 

x* n -x n = (y + E<p(x* n )/a n 
x n+ i -x* n = (y + e(f>(x* n )/a n+1 , (9) 

where <p(x) is the corrugation function introduced for the equation (1) with an amplitude 
of modulation e <C yo- We take into account the small-angle limit of ([7]) tana a. The 
corresponding relation between a n and a n+ i is 

a n+ i = a n -2(3. (10) 

The angle (3 is defined at the point x* n by a tangent line for <p(x) such that tan(/?) = 
Ed(f)(x) / dx\ x=x » (see Fig. lb). Redefining z n = x* n and neglecting E(fi(x*)/y = e(p(z)/y due 
to the conditions e/yo -C 1 and \4>(z)\ < 1, we obtain from (Q that the new relation reads 

z n +i = x n +i + yo/a n+ i = z n + 2y /a n+ i. (11) 
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Eventually, the following Hamiltonian map describes ray dynamics 

a n+ i = a n - 2e<p'(z n ) 

z n +i = z n + 2y /a n+1 , (12) 

where <j> r (z) = d(p(z)/dz. It is simple to see that \d(a n+ i, z n+ i)/d(a n , z n )\ = 1. In the case 
when d>'(z) is periodic, not necessarily differentiable, the map ffl~2l) is called the Ulam map 



20| | . This map corresponds to a ray propagation when any two nearest bouncings occur 
from the opposite boundaries. The necessary condition ensures this effect when the maximal 
angle determined the corrugation k^s/ft is smaller than the minimum momentum a m i n = h. 
Hence, the condition of the validity of the map (iTTj) is 

2e/\ < 7c/k y . (13) 

In what follows, we consider 

4>(z) = cos(k z). 

For chaotic dynamics the angles are equally distributed after some time in some chaotic 
region of the phase space for any single initial condition chosen in the chaotic domain (see a 
phase portrait in Fig. 2). The phase portrait is obtained for a single trajectory by iteration 
of the map (fl2"|) when the following variable change is made: u = a/2koE and ip = k$z. In 
this case, the chaos control parameter is 

M = y Q /e, (14) 



and it determines the maximal size of the chaotic region \u\ < yM. For the angles, the 
maximal size of the chaotic region is a max ~ 2k y yje/y$. 



7 



A CANONICAL VARIABLE CHANGE 



Quantum or wave interference leads to localization of chaos and, consequently, to the 
essential difference of the initial profile spreading from the ray dynamics 



2l|. A quantum 



mechanical counterpart o: 
acceleration mechanism 



the map f fT2l was the subject of many studies, related to the Fermi 



19 



22 



231 ]. Some examples of the quantum Fermi acceleration 



dynamics have been considered in the framework of the so-called Generalized Canonical 



Transformation 



221. It shou 



time" canonical variables 



d be admitted that the Ulam map is written in the "energy- 



201 ]. including the map (1T21 . Since the longitudinal coordinates 



play a role of time, we introduce the time parameter in the following dimensionless form 
vt = koz, where v = koyo, while the energy is a, and it scales by 2: a/2 — > a. Therefore the 
Hamiltonian equations for map (fT2"|) in the new variables are 



da/ dy = k e sin vt Yln=-oo $ (v ~ n 
dt/dy = 1/a, 



(15) 



where the formal time parameter y is the dimensionless transversal coordinate, which is 
scaled by 2y . The Hamiltonian is 



if = ln|a| cosvt 5{y — n). 



Vo 



(16) 



Quantization of the system (jlBj) is convenient to carry out in the framework of the transver 



24 



261 ] . Moreover, this quantization 



sal momentum and coordinate (p, y) canonical pair 
is natural. In the limit a < 1 a simple variable change could be suggested. Namely, the 
Hamiltonian ffTBT) can be rewritten in a new form, such that t will be the real time pa- 



rameter, while the dimensionless transversal coordinate y and the corresponding transversal 
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momentum p will be canonical variables 25j . Inverting the second equation in (|T5|) we obtain 

y = a. (17) 
Now, let us introduce the new unperturbed Hamiltonian TCo(p), such that 

dHo(p)/dp = a, da/dHo = a. (18) 

It follows from that 

da/dp = (da/dHo) ■ (dH (p)/dp) = a 2 . (19) 

Solving fllBJ and we obtain that 

a = —l/p and H-oip) = In \a\ = In \l/p\. (20) 

This transformation for the unperturbed system corresponds to the mapping of the interval 
a G (0,a max ) to the interval p G (— oo, — l/a max ) and the interval a G (— a ma x,fy to 
p G (l/a>max, oo). From the first equation in f|T5l) we have 

ryn+i-0 ^ a rVn+i-0 

/ -j-dy = k e / sin ut5(y — n)dy. (21) 

The integrands from both sides of (l2Til can be transformed. Using (|T71) and ( 1T91) we have for 

the left-side integrand the following expression 

da dt da . &Ho . , . 

dy dy dp ^ dp ^ 



By use ( ITT]) , the right hand side integral (rhsi) can be rewritten approximately in the form 

'Vn+l—O pt n -{-i—0 
'Vn-0 

Then taking into account that for the almost equidistant part o 
we have &Ho(p)/dp ~ const, and one can write approximately 27| that 



rVn+i—O rtn+i—0 

/ sin vt5{y — n)dy = I sin vtb it — t n )dt. (23) 

Jyn-0 Jt n —0 



the spectrum, when p 1, 



^(t-t n )^2/ 2 ^%-n). (24) 



Hence the rhsi in (12 ip and (123|) reads 

PVn+l-O J rtn+l-0 ^ 

/ sin z/t5(?/ — n)dy m — / cos vt— 5 (t — t n )dt 

Jy„-0 v Jtn-0 dt 

i d 

« — / y cos ut— S(y -n)dy. (25) 

Then we obtain from (T2TT ) — f|24T) the following equations of motion 

k e °° d 

p = cosz4 } —5(y — n) an d y — —l/p (26) 

^-^ ay 

n=—oo 

with the Hamiltonian 

, oo 

Ti = - In |p| + — cos z4 ^ ~ n ) = ^o(p) + t). (27) 

n=— oo 

QUANTIZATION OF THE ULAM MAP 

A quantization procedure now corresponds to the consideration of the Shrodinger equa- 
tion for the wave function. The semiclassical consideration requires that the width of the 
perturbative potential V in (12?]) is larger than the wavelength. Therefore, it is necessary to 
restrict the summation in the Fourier expansion of the <52tt(#) potential. Thus we obtain a 
new potential with the width of a spike equaled to 2n/N: 

N N 



5 N (9) = 2^2 cos k9 -1 = expikO, 



k=0 k=-N 

where 5n(0) tends to £271- (#) at iV tends to infinity. 
Floquet theory 

Since the Hamiltonian (12 7p is periodic in time, the Floquet theory is used in what follows. 
In this case the wave function due to the Floquet theorem is 

i;(y,t) = e- lAt My,t), 
10 



where if)\(y, t + 2ir/v) = ^\(y, t) is the periodic eigenfunction with the period of the pertur- 
bation. The Schrodinger equation 

ih^(y,t)=H(p,y,i)Tl>(y,t) (28) 

corresponds to the eigenvalue problem for the Floquet operator F 

i^A = hAip A . (29) 

The Floquet operator is 

P = ~ih^ + K(p, y, t) = F + V(y, t). (30) 

The solution of (J29~j) is considered as a superposition of the unperturbed basis 

\nj) = \n)\j) = ^e a ^e-» vt , (31) 

which is the eigenfunction of the unperturbed system F Q = —ihd/dt + 'H andp|n) = hn\n). 
Taking into account that the perturbation is periodic in y, we obtain for the eigenfunction 

Mv,t) = \Mv,t)) = e' Wy J2Mn,f)\n,j), (32) 

where 9 is a quasi-momentum. The coefficients of the expansion <f> n j = <p\ } e(n,j) = 
(j,n\ipA(y,t)) can be found from the following equation 

S(n + d)- hvj 

n' 

where £{n + 9) = (n\Ti (p + 6)\n) and , <j) n i^ < oo. 
An exact solution in the linear approximation 

Some simplification of the equation can be made. Taking into account that in the range 
of the almost equidistant spectrum for p 1, one can consider approximately that £(po + 
9 + hn) ^£ + hv(p )n, where £ = £(p + 9) and v(p ) = ^-£(po + d). 
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K e 
2v 



} n',j+l 



+ 



>n',j-l\ 



hA(j) n . 



3i 



(33) 



To lift the summation over n' in (1331). a sine -function sincx = is introduced. It 



possesses 5-like properties 



2_^sincir(m — a) — 1, and sincO = 1. 

m 

Then the solution of 0331) is cast in the form of the sinc-function 



(34) 



b n>j = Zj8inc7r[n - (A - £ )/v\, 



(35) 



where coefficients Zj correspond to the following equation 



hvn — hvj — h(A — So/h) Zjsmc ir[n — (A — £q)/v 



+^ [Z j+1 + = 0. 



Counting for the fixed n and j that the quasienergy spectrum is 



hA = £ + /ira;, 



we obtain the following relation for the Bessel functions 28] 



2jZ j =^-[Z j+1 + Z j . 1 ] 
hv 



(36) 



(37) 



(38) 



with the solution Zj = Jj{¥^) = Jj{^)- m the case when k y 3> 1, the effective number of 
transitions due to the perturbation is restricted by the relation for the Bessel function of a 



small argument 2e/\<^ 1 



such that 



^(f)~(f) J =exp{-jln(A/2,)}. 



Therefore, the effective number of interacting eigenmodes is restricted by the relation 



(39) 



\j\ < j ~l/ln(A/2e). 



(40) 
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In the opposite case when the period of the corrugation is larger than the transversal size 
of the channel, one obtains from (Tl3|) that 2e/A ^> 1. In this case the effective number of 
transitions due to the perturbation is restricted by the relation 

\j\ <j ~2e/\, (41) 



28] 



since the Bessel functions in this case decay faster than the exponential when j > j 
One should also recognize that in a general case the approximation for the quasi-equidistant 
spectrum could be not valid. Then the approximate analysis should be carried out beyond 



the linear approximation used here for the expression (|4"II) 25]. 



CONCLUSION 

The ratio between the wavelength and the amplitude of the corrugation is important not 
only for the validity of the Ulam map (fl~2l . but because it also determines the relation on the 
ray chaos for the quantum or wave process. Since chaotic ray dynamics is strongly damped 
by quantum effects, the quantum system behaves like its classical counterpart of eq. (TT2"|) 
on a finite time scale only. Therefore, one can discuss the chaos localization phenomenon 
by quantum effects. It means that quantum chaotic dynamics is localized in some energy 
scales determined by ff4T3l) and f|4Tj) . Borrowing some familiar terminology, one defines this 
range as a localization length. In the units of the unperturbed spectrum, the localization 
length is An < j^v/v. For the angles it reads 

1 i {k e) ■ a : A < e 

Aa q = = « { ( 42 ) 

p ° (p ° +/iAn) ((M-i-N^r-o : A» £) 

and the angles' width Aa q also defines the width of the spreading wave packet. We cannot 
but to admit two important characteristics (or parameters) which determine the localization 

13 



lengths in ( 1421) . The first one follows from (1B1) . (TT41 and ( 14T1) . The semiclassical parameter h, 
the localization length 2e/X and the chaos control parameter M form the following relation, 
familiar in quantum chaos [17] 

l/M=(A/23/ ).(2e/A). (43) 

The conditions of quantum chaos are X/yo <C 1 and M 3> 1. This interplay of classical 
or ray chaos with the quantum interference is an important mechanism of quantum chaos 



or quantum localization of classica 
Rayleigh condition k e < 0.38 



chaos. The second one is concerned with the so-called 
29j imposed on the system for a spreading wave. In our 
case koE < h -C 1, which follows from (|T3|) . One may also understand this result from a 
point of view of wave dynamics. The excitation of an initial wave with an arbitrary angle 
ao means that an initial wave packet is a superposition of eigenmodes of a m (see (jZJ)). If 
the corrugation is weak enough, the wave packet is a narrow superposition of the Floquet 
eigenfunctions of (j3J) with a m being close to the initial angle a , and its width Aa q remaining 
a constant value during the evolution. An estimation of this width is determined by Eq. 
( 1421) : \a m — ao\ < Aa q . We should stress that the present analysis can be considered as a 
qualitative estimation of this width as well. 
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Figure captions 

Fig. 1. A ray dynamics sketch for constructing map ffTUl) : a' = a — 2(3. 
Fig. 2. Phase portrait of the map (fT2"j) for a single trajectory with an initial condition 
(if),u) = (0, 0.1). The following variable change is made u = ajlk^e and if) = k^z, while the 
chaos control parameter is M = yo/e. The trajectory is shown after 125000 iterations for 
M = 2vr • 125. 
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